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PHYSIOLOGICAL PROFILING 

Cross Reference To Related Applications 

This application claims priority from U.S. Provisional Application Serial No. 
60/234 023, filed on September 20, 2000. 

Statement as to Federally Sponsored Research 
Funding for the work described herein was provided by the federal government, 
which may have certain rights in the invention. This work was supported by National 
Heart, Lung, and Blood Institute Grant 1P50-HL-54998. 

BACKGROUND 

L Technical Field 

The invention relates to methods and materials involved in identifying 
relationships among physiological determinants (parameters) associated with complex 
physiological processes that contribute to normal and pathological states of an organism. 

2. Background Information 

Genetic studies of complex multifactorial diseases such as asthma, hypertension, 
non-insulin-dependent diabetes mellitus (NIDDM), and insulin-dependent diabetes 
mellitus (IDDM) remain challenging due to heterogeneity in the clinical presentation of 
these diseases among patient populations. In addition, the modest contribution of each 
gene and/or the study of phenotypes that are distant from these gene effects, or both, have 
made identifying genes involved in these diseases difficult. Difficulties in elucidating the 
genetic basis of multifactorial diseases have become apparent from results obtained from 
total genome scans for quantitative trait loci (QTL) associated with asthma, hypertension, 
NIDDM, and IDDM in diverse human populations. (See Bleecker et al (1997) Am J 
RespirCriL Care Med. 156:S1 13-6; Julier et al (1997) Hum Mol Genet 6:2077-85; 
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Kmshkal et al (1999) Circulation 99:1407-1410; Hanis etal (1996) Nat Genet 13:161- 
166; and J. A. Todd (1995) Proc Natl Acad Sci USA 92, 8560-8565). 

Furthermore, although genome-wide scans directed at the genetic basis of 
hypertension in rats have identified rough locations of genes on almost every rat 
chromosome, with loci confirmed on chromosomes 1, 2, 3, 5, 10, and 13 (J. P. Rapp 
(2000) Physiol. Rev. 80:135-172), no actual genes have been identified. The need for 
improved analytical tools, in addition to better phenotyping protocols, for identifying 
genes influencing complex phenotypes has been well articulated by Nadeau and Frankel 
(Nadeauefa/. (2000) Nat. Genet. 25,381-384). 

SUMMARY 

The invention provides methods and materials related to identifying relationships 

□ 

y3 among physiological traits - herein referred to as "physiological determinants." More 

^ specifically, the invention provides a new analytical procedure for identifying 

£3 1 5 relationships among physiological determinants associated with complex physiological 

fjj 

y processes that contribute to normal and pathological states of an organism. The 

analytical procedure, termed "physiological profiling," involves, in broad form, three 
£3 steps. First, a set of physiological determinants is identified. Second, correlation values 

p j are determined between pairs of physiological determinants for all possible pairs within 

l 20 the set. Third, the correlation values are organized into a clustered correlation matrix by 
organizing the corresponding physiological determinants along the axes of the matrix 
using a clustering method. From the resulting "physiological profile," relationships 
between determinants can be identified. Physiological profiling can be used to 
characterize physiological processes in normal and diseased organisms. Results of 
25 physiological profiling can be used to classify diseased and/or normal organisms into 

groups based on correlation patterns determined. Physiological profiling also can be used 
in conjunction with genetic linkage analysis or gene expression profiling for functional 
genomics studies or clinical diagnosis. 

In one embodiment, the invention provides a method of identifying relationships 
30 among physiological determinants within a set of physiological determinants. The 
method involves (1) determining a correlation value between two physiological 
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determinants for all possible pairs of physiological determinants within the set; (2) 
constructing a correlation matrix using the determined correlation values; (3) constructing 
a clustered correlation matrix from the correlation matrix by clustering physiological 
determinants using a clustering method, and (4) identifying relationships among 
5 physiological determinants from the clustered correlation matrix. The clustering method 
can be based on known physiological relationships, known genetic linkages, or gene 
expression profiles. Alternatively, the clustering method can be a statistical method that 
does not rely on known physiological relationships, genetic linkages, or gene expression 
profiles. 

1 0 In another embodiment, the method can involve constructing a colored clustered 

correlation matrix using a plurality of colors such that each color indicates a selected 
degree of correlation. The patterns of colors in the clustered correlation matrix can be 
used to identify physiological relationships. 

In another embodiment, the set of physiological determinants can include at least 

p 15 10, 20, or 50 determinants. 

fU 

y 
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In another embodiment, the first member of each pair of physiological 
determinants can be derived from an individual and the second member of each pair of 
physiological determinants is the mean of physiological determinants from a population 
of individuals; and the correlation value is determined by a method that includes 



yj 

□ 20 measuring the difference between the first member and the second member. 
0 



In another embodiment, the invention provides a method of assessing the 
physiological response of an organism to a challenge. The method includes a first, 
second, and third step. The first step involves constructing a first clustered correlation 
matrix using a set of physiological determinants. The first set of correlation values for all 
25 pairs of determinants in the set is obtained prior to the challenge. The second step 
involves constructing a second clustered correlation matrix using the same set of 
physiological determinants, and the second correlation values for all pairs of determinants 
in the set are obtained during or subsequent to the challenge. The third step involves 
comparing the first and second clustered correlation matrices to assess the physiological 
30 response of the organism to the challenge. The challenge can be, for example, a drug 
administration, an allelic substitution, or an environmental stressor. 
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In one embodiment, the correlation values in the matrices are represented by a 
plurality of colors, each color indicating a selected degree of correlation. In another 
embodiment, multiple clustered correlation matrices can be compared by comparing the 
patterns of colors of each matrix. 
5 In another embodiment, the invention provides a method of assessing the change 

in physiological state of an organism or organisms over time. This method includes a 
first, second, and third step. The first step involves constructing a first clustered 
correlation matrix using a set of physiological determinants. The correlation values for 
all pairs of determinants in the set are obtained at a first time point. The second step 

10 involves constructing a second clustered correlation matrix using the same set of 

physiological determinants. The correlation values for all pairs of determinants in the 
second step are obtained at a second time point. The third step is comparing the first and 
second clustered correlation matrices to assess the change in physiological state of the 
organism from the first to the second time point. In another embodiment, more than two 

1 5 time points can be compared in this manner. The correlation values can be represented 
by a plurality of colors with each color indicating a selected degree of correlation. The 
clustered correlation matrices are compared by comparing the patterns of colors in the 
clustered correlation matrices. 

In another embodiment, the invention provides a method of partitioning 

20 organisms into homogeneic subclasses. The method involves comparing the 
physiological profiles of the organisms and then partitioning the organisms into 
homogeneic subclasses based on differences in the physiological profiles. In one 
embodiment, expression profiling can be used to further partition the organisms into 
additional homogeneic subclasses based on expression profiling results. In another 

25 embodiment, the organisms can exhibit a multifactorial disease condition. 

In another embodiment, the invention provides a method of assigning an organism 
to a homogeneic subclass of organisms. The method includes generating a physiological 
profile of the organism and identifying the organism as belonging to a homogeneic 
subclass based on the physiological profile. In another embodiment, the homogeneic 

30 subclass of organisms exhibits a multifactorial disease condition. 

In another embodiment, the invention provides a method of determining the 
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contribution(s) of a gene or genes to a physiological process in an organism. The method 
involves a first, second, third, and fourth step. The first step involves generating a first 
expression profile and a first physiological profile of the organism before a challenge. 
The second step involves generating a second expression profile and a second 
5 physiological profile of the organism during or after the challenge. The third step 
involves comparing the first expression profile and first physiological profile with the 
second expression profile and second physiological profile. The gene or genes are 
identified by the difference or differences in the first and second expression profiles. The 
physiological contributions of the same gene or genes are indicated by changes in the first 
1 0 and second physiological profiles. 

In another embodiment, the invention provides a method of determining the 

□ contribution(s) of a gene or genes to a physiological process in an organism. The method 
% involves a first, second and third step. The first step is generating a first expression 

EH profile and a first physiological profile of the organism at a first time. The second step is 

□ 

n j 1 5 generating a second expression profile and a second physiological profile of the organism 
*** at a second time. The third step is comparing the first expression profile and first 

s physiological profile with the second expression profile and second physiological profile. 

£3 

The gene or genes are identified by the difference or differences in the first and second 
jj; expression profiles. The physiological contributions of the gene or genes are indicated by 

□ 20 changes in the first and second physiological profiles. 

r In another embodiment, the invention provides a computer-readable medium that 

includes a physiological profile. In another embodiment, the physiological profile has a 
plurality of colors, each color indicating a selected degree of correlation. 

In another embodiment, the computer-readable medium can have stored 
25 computer-readable instructions for performing the above-described methods. 

In another embodiment, the invention provides a method of determining whether 
a hypertensive patient is a modulator or non-modulator. The method involves 
determining the allelic status of a gene encoding renin in the patient. The allelic status of 
a patient is determined by identifying which allele of a relevant gene, among various 
30 possible alleles of that gene, is possessed by that patient. 

In another embodiment, the invention provides a method of determining whether 
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a patient is at risk for hypotension following administration of a vasoconstrictor agent. 
The method involves determining the allelic status of a gene encoding NOSII in the 
patient. 

In another embodiment, the invention provides a method of determining whether 
a patient is at risk for hypotension following administration of a vasoconstrictor agent. 
The method involves determining the allelic status of a gene encoding NOSIII in the 
patient. 

In another embodiment, the invention provides a method for modifying or 
supplementing actuarial tables for life and health insurance. The method involves 
identifying homogeneic subclasses of organisms, e.g. humans, as described earlier, and 
modifying or supplementing actuarial tables based on the identified homogeneic 
subclasses. 

Unless otherwise defined, all technical and scientific terms used herein have the 
same meaning as commonly understood by one of ordinary skill in the art to which this 
invention pertains. Although methods and materials similar or equivalent to those 
described herein can be used in the practice or testing of the present invention, suitable 
methods and materials are described below. All publications, patent applications, 
patents, and other references mentioned herein are incorporated by reference in their 
entirety. In case of conflict, the present specification, including definitions, will control. 
In addition, the materials, methods, and examples are illustrative only and not intended to 
be limiting. 

Other features and advantages of the invention will be apparent from the 
following detailed description, and from the claims. 

DESCRIPTION OF DRAWINGS 

Figure 1 is a comprehensive linkage map of 81 determinant phenotypes (96 QTL) 
in the autosomal genome of F2 male progeny (n=l 13) from an SS/JrHsd/Mcw and 
BN/SsNHsd/Mcw intercross. Vertical bars on the left side represent the 95 % confidence 
intervals (CI) of individual QTL. Green bars indicate CI from parametric analysis, while 
orange bars indicate CI from non-parametricWalysis. Phenotype designations and peak 
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LOD scores (green = parametric) and Z-scores (orange = non-parametric), respectively, 
are presented on the right of eaoji chromosome. 

Figure 2 is a randomized colored correlation matrix of BN phenotypes. Strong 
positive correlations are represented in red, strong negative correlations are blue, while 
5 low correlations are in gray and black. 

Figure 3 is a physiological profile of BN phenotypes ordered by functional 
clustering using Guyton's model of blood pressure control. Strong positive correlations 
are represented in red, strong negative correlations are blue, while low correlations are in 
gray and black. 

10 Figure 4 is a composite matrix of two physiological profiles, one generated using 

functional clustering and the second generated using purely statistical clustering. Strong 
positive correlations are represented in red, strong negative correlations are blue, while 
j3 low correlations are in gray and black. 

Figure 5 is two physiological profiles consisting of phenotypes associated with 
C3 1 5 regulation of blood flow for parental BN and SS rats. Strong positive correlations are 
\ |1 represented in red, strong negative correlations are blue, while low correlations are in 

H£ gray and black. 

P Figure 6 is two composite physiological profiles of parental BN and F2 progeny 

rats generated by overlaying functionally clustered correlation matrices with algorithm 
£3 20 clustered correlation matrices. 

II Figure 7 A is a comparison of the physiological profile of all F2 progeny rats with 

the physiological profile of progeny rats that fall in the left 10 % tail of a distribution 
after a salt challenge. Strong positive correlations are represented in red, strong negative 
correlations are blue, while low correlations are in gray and black. 

25 Figure 7B is a comparison of the physiological profile of all F2 progeny rats with 

the physiological profile of progeny rats that fall in the right 10 % tail of a distribution 
after a salt challenge. Strong positive correlations are represented in red, strong negative 
correlations are blue, while low correlations are in gray and black. 

Figure 8 A is a physiological profile of phenotypes associated with arterial blood 

30 pressure in F2 male rats homozygous SS for DIOMghH (NOSII gene region). Strong 
positive correlations are represented in red, strong negative correlations are blue, while 
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low correlations are in gray and black. The expanded insert represents correlations 
among blood pressures determined immediately before, during, and after administration 
of norepinephrine, angiotensin II, and acetylcholine. 

Figure 8B is a physiological profile consisting of phenotypes associated with 
5 arterial blood pressure in F2 male rats homozygous BN for D10Mghl4 (NOSH gene 
region). Strong positive correlations are represented in red, strong negative correlations 
are blue, while low correlations are in gray and black. The expanded insert represents 
correlations among blood pressures determined immediately before, during, and after 
administration of norepinephrine, angiotensin II, and acetylcholine. 
10 Figure 9A is a graph illustrating the correlation between mean arterial pressure 

before and after infusions of norepinephrine in F2 rats homozygous SS (open circles) for 
the NOSII gene and those homozygous BN (closed circles) for the NOSII gene. 

Figure 9B is a bar graph summarizing the average levels of mean arterial pressure 
before (solid bars) and following completion (open bars) of the intravenous infusions of 
1 5 three doses of norepinephrine in male rats carrying the SS or BN allele at NOSII. 

Figure 10 is two physiological profiles of French Canadian and African American 
hypertensive patients. Strong positive correlations are represented in red, strong negative 
correlations are blue, while low correlations are in gray and black. 

20 DETAILED DESCRIPTION 

The invention provides methods and materials related to identifying relationships 
among physiological traits - herein referred to as "physiological determinants." More 
specifically, the invention provides a new analytical procedure for identifying 
relationships among physiological determinants associated with complex physiological 

25 processes that contribute to normal and pathological states of an organism. The 

analytical procedure, termed "physiological profiling," involves, in broad form, three 
steps. First, a set of physiological determinants is identified. Second, correlation values 
are determined between pairs of physiological determinants for all possible pairs within 
the set. Third, the correlation values are organized into a clustered correlation matrix by 

30 organizing the corresponding physiological determinants along the axes, for example, 
top, bottom, or sides of the matrix using a clustering method. From the resulting 
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"physiological profile," relationships between determinants can be identified. As used 
herein, the term "physiological profile" refers to a clustered correlation matrix generated 
using (1) a set of physiological determinants, ordered using a clustering method, and (2) 
the correlation values determined for all possible pairs of physiological determinants in 
5 the set. As used herein, the term "physiological profiling" refers to an analytical 

procedure involving (1) identifying a set of physiological determinants, (2) determinating 
correlation values for all possible pairs of determinants with the set, and (3) generating a 
clustered correlation matrix by organizing the correlation values into a matrix using a 
clustering method that orders the determinants in a non-random fashion along the axes of 
10 a matrix. Physiological profiling can be used to characterize physiological processes in 
normal and diseased organisms. Results of physiological profiling can be used to classify 
P diseased and/or normal organisms into groups based on correlation patterns determined. 

Physiological profiling also can be used in conjunction with genetic linkage analysis or 
gene expression profiling for functional genomics studies or clinical diagnosis. 

S 15 

'% Physiological determinants 

» b The first step in generating a physiological profile is identification of a set of 

J3 physiological determinants. As used herein, the term "physiological determinants" refers 

% to physiological traits that can be determined experimentally or derived from 

C3 20 experimentally measured data. For example, a measured physiological determinant can 
be weight (e.g. weight of an organism or an organ such as a kidney), volume (e.g. urine 
volume), or blood pressure (e.g. diastolic or systolic blood pressure). A derived 
physiological determinant can be, for example, mean blood pressure, standard deviation 
of mean arterial pressure, or the difference between (1) blood flow/gram of kidney weight 
25 after administration of a drug and (2) control blood flow per gram of kidney weight. 

Physiological determinants can be obtained before, during, or after a challenge. 
A challenge can be any condition or event that triggers a physiological response or alters 
homeostasis. The challenge can be, for example: a disease condition, one or more allelic 
substitutions, an environmental stressor (e.g. hypoxia, high salt intake), contact with a 
30 naturally- or non-naturally occurring chemical or macromolecule, infection by a 
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biological material (e.g. bacteria, viruses, prions), and the presence or absence of 
exercise. 

One example of a derived physiological determinant is "delta renal blood flow 
from Angiotensin II dose 2 minus control renal blood flow." In this example, the 
5 physiological determinant is derived by subtracting renal blood flow determined before a 
challenge, from delta renal blood flow determined after a challenge. The challenge is 
Angiotensin II. 

Physiological determinants reflect the status of the relevant complex 
physiological system, for which the determinants serve as estimates of biological 

1 0 function. Complex physiological systems include, without limitation, the respiratory 
system, cardiovascular system, nervous system, digestive system, endocrine system, 
immune system, lymphatic system, renal system, skeletal system, catabolic and metabolic 
systems, and the digestive system. 

Physiological determinants can include coronary determinants associated with 

1 5 mechanical, electrical, and biochemical functions in the heart, and with the heart's ability 
to resist ischemia. Examples include, without limitation, ischemic peak contracture 
(mmHg), ischemic time to onset of contracture (sec), ischemic time to peak contracture 
(sec), post-ischemic coronary flow rate (mL/min), enzyme leakage (IU/g wet weight), 
heart rate (beats/min), infarct size (% LV), left ventricle developed pressure (mmHg), left 

20 ventricle diastolic pressure (mmHg), left ventricle systolic pressure (mmHg), recovery 
coronary flow rate (% recovery), recovery developed pressure (% recovery), recovery 
heart rate (% recovery), recovery systolic pressure (% recovery), coronary flow rate 
(ml/min/g), enzyme leakage (IU/g wet weight), post-ischemic heart rate (beats/min), pre- 
ischemic heart wet weight (g), left ventricle developed pressure (mmHg), left ventricle 

25 diastolic pressure (mmHg), and pre-ischemic left ventricle systolic pressure (mmHg). 

Physiological determinants can be associated with the vascular system and 
include vascular responsiveness to acute vasoconstrictors and dilators, vascular function, 
and the susceptibility to developing injury in response to a high salt diet. Examples 
include, without limitation, dilator response to acety choline EC50 (1 x 10~ 7 mole), dilator 

30 response to acety choline Log EC50 (Log molar), fast slope of phenylephrine-induced 
contraction (gram/min), maximum force (g) per wet weight of aorta (gram/min), % 
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maximum relaxation acetylcholine (%), % maximum relaxation of phenylephrine- 
induced contraction by 0 % O2 (%), % maximum relaxation of phenylephrine-induced 
contraction by 10 % O2 (%), % maximum relaxation of phenylephrine-induced 
contraction by 5% O2 (%)> % maximim relaxation sodium Nitroprusside (%), constrictor 
response to phenylephrine EC50 (1 x 10' 7 mole), constrictor response to phenylephrine 
Log EC50 (Log molar), dilator response to sodium nitroprusside EC50 (1 x 10" 7 mole), 
dilator response to sodium nitroprusside Log EC50 (Log molar), and slow slope of 
phenylephrine-induced contraction (gram/min). 

Physiological determinants can be associated with renal function such as blood 
pressure responsiveness to acute vasoconstrictors and dilators, and renal tubular function 
and susceptibility to developing renal injury in response to high salt diets. Examples 
include, without limitation, baseline HR for Angll dose-response relationship 
(beats/min), NE dose-response relationship (beats/min), baseline MAP for Angll dose- 
response relationship (mmHg), and baseline MAP for NE dose-response relationship 
(mmHg). Examples also include high salt creatinine clearance (mL/min), low salt 
creatinine clearance (mL/min), delta HR to 10 ng/kg/min Angll (beats/min), delta HR to 
0.2 ug/kg/min NE (beats/min), delta HR to 25 ng/kg/min Angll (beats/min), delta HR to 
0.5 ug/kg/min NE (beats/min), delta HR to 50 ng/kg/min Angll (beats/min), delta HR to 
1 .0 ug/kg/min NE (beats/min), delta HR to 5 ng/kg/min Angll (beats/min), and delta HR 
to 0.1 ug/kg/min NE (beats/min). Examples also include change in heart rate with salt 
depletion (beats/min), high salt heart rate (beats/min), low salt heart rate (beats/min), pre- 
to post-control delta HR following ANGII (beats/min), pre to post control delta HR 
following NE (beats/min), delta MAP to 10 ng/kg/min Angll (mmHg), delta MAP to 0.2 
ug/kg/min NE (mmHg), delta MAP to 25 ng/kg/min Angll (mmHg), delta MAP to 0.5 
ug/kg/min NE (mmHg), delta MAP to 50 ng/kg/min Angll (mmHg), delta MAP to 1.0 
ug/kg/min NE (mmHg), delta MAP to 5 ng/kg/min Angll (mmHg), delta MAP to 0.1 
ug/kg/min NE (mmHg), change in mean arterial pressure with salt depletion (mmHg), 
high salt mean arterial pressure (mean of three days of high salt pressure recordings) 
(mmHg), low salt mean arterial pressure (one day of recording following salt depletion) 
(mmHg), pre- to post-control delta MAP following ANGII (mmHg), pre- to post-control 
delta MAP following NE (mmHg), high salt plasma creatinine (mg/dL), low salt plasma 
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creatinine (mg/dL), change in plasma renin activity with salt depletion (ls-hs) (ng 
angl/mL/hr), high salt plasma renin activity (ng angl/mL/hr), low salt plasma renin 
activity (ng/mL/hr), high salt urinary excretion of sodium (mEq/day), low salt urinary 
excretion of sodium (mEq/day), high salt urine microalbumin excretion (mg/day), high 
5 salt urine osmolality (mOsm/L), low salt urine osmolality (mOsm/L), and high salt urine 
protein excretion (mg/day). 

Physiological determinants also can be associated with lung functions such as 
airway methacholine sensitivity, pulmonary vascular mechanics, pulmonary endothelial 
angiotensin converting enzyme activity, and pulmonary endothelial redox status in 
10 normal and chronically hypoxic conditions. Examples include, without limitation, alpha, 
(a statistical measure of characterizing the white noise component of blood pressure - 
1/mmHg), body weight (kg), FAPGG metabolism-surface area product (mL/min x kg), 
ji lung dry weight/body weight ratio (g/kg), hematocrit (%), MB+ metabolism-surface area 

product 1 (mL/min x kg), MB+ MSAP 3 (mL/min x kg)/FAPGG MSAP (mL/min x kg), 
C5 1 5 MB+ metabolism-surface area product 3 (mL/min x kg), methacholine ED50 (mg/kg), 
rl right ventricle/Left ventricle weight ratio (w/w ratio), and r @flow = lOOmL/min/g ("r" 

represents left ventricular resistance, mmHg x min x kg/ml). 
q Physiological determinants can be associated with respiration such as respiratory 

Si control mechanisms and the pattern of breathing and lung function in the conscious state 

□ 20 under acute conditions of hypoxia, hypercapnia, and exercise. Examples include, without 
rl limitation, heart rate during control (co) HYPERCAPNIA (beats/min), heart rate during 

control (co) HYPOXIA (beats/min), change in heart rate from rest to run (delta re v rn) 
(beats/min), change in heart rate from rest to walk (delta re v wk) (beats/min), heart rate 
during minute 7 of hypercapnia (b2) (beats/min), heart rate during minute 7 of hypoxia 
25 (b2) (beats/min), heart rate treadmill resting 3 minute average (re) (beats/min), heart rate 
running 30 second average (rn) (beats/min), and heart rate walking 30 second average 
(wk) (beats/min). Examples also include mean arterial pressure during control (bl) 
HYPERCAPNIA (mmHg), mean arterial pressure during control (bl) HYPOXIA 
(mmHg), change in mean arterial pressure from rest to run (delta re v rn) (mmHg), 
30 change in mean arterial pressure from rest to walk (delta re v wk) (mmHg), mean arterial 
pressure during minute 7 of hypercapnia (b2) (mmHg), mean arterial pressure during 



12 



Docket 13482-0020 

minute 7 of hypoxia (b2) (mmHg), mean arterial pressure treadmill resting 3 minute 
average (re) (mmHg), mean arterial blood pressure treadmill 30 second average (rn) 
(mmHg), and mean arterial blood pressure treadmill 30 second average (wk) (mmHg). 
Examples also include control (co) PaC0 2 HYPOXIA (mmHg), change in PaC0 2 
5 between Control (co) and hypoxia (h2) (mmHg), hypoxia (h2) PaC02 (mmHg), arterial 
PCO2 at rest 30 second average (re) (mmHg), arterial PCO2, running 30 second average 
(rn) (mmHg), arterial PCO2 walking 30 second average (wk) (mmHg), control (co) PaC>2 
HYPOXIA (mmHg), change in Pa02 between Control (co) and hypoxia (h2) (mmHg), 
hypoxia (h2) Pa02 (mmHg), arterial PO2 at rest 30 second average (re) (mmHg), arterial 

10 P0 2 running 30 second average (rn) (mmHg), arterial P0 2 walking 30 second average 
(wk) (mmHg), change in PCO2 from rest to run (delta re v rn) (mmHg), and change in 
PCO2 from rest to walk (delta re v wk) (mmHg). Examples also include control (co) pH 
HYPOXIA (pH), change in pH between Control (co) and hypoxia (h2) (pH), change in 
pH from rest to run (delta re v rn) (pH), change in pH from rest to walk (delta re v wk) 

1 5 (pH), hypoxia (h2) pH (pH), arterial pH at rest 30 second average (re) (pH), arterial pH 
running 30 second average (rn) (pH), arterial pH walking 30 second average (wk) (pH), 
change in PO2 from rest to run (delta re v rn) (mmHg), and change in PO2 from rest to 
walk (delta re v wk) (mmHg). Examples also include rectal temperature for control 
HYPERCAPNIA (°C), rectal temperature for control HYPOXIA (°C), change in rectal 

20 temperature from rest to post exercise (°C), change between control and hypercapnic 
rectal temperature (°C), change between control and hypoxic rectal temperature (°C), 
rectal temperature following running on treadmill (°C), rectal temperature after 
hypercapnia (°C), rectal temperature after hypoxia (°C), rectal temperature at rest (°C), 
pulmonary ventilation (VE) control (co). HYPERCAPNIA (mL/min), and pulmonary 

25 ventilation (VE) control (co). HYPOXIA (mL/min). Examples also include % change 
from (co)_in ventilation to (h2) hypercapnia, % change from (co) in ventilation to (h2) 
hypoxia, pulmonary ventilation (VE) at hypercapnia from minute 2-3 (hi) (mL/min), 
pulmonary ventilation (VE) at hypoxia from minute 2-3 (hi) (mL/min), % change from 
(co)_in ventilation to (h2) hypercapnia, % change from (co) in ventilation to (h2) 

30 hypoxia, pulmonary ventilation (VE) at hypercapnia from minute 9-10 (h2) (mL/min), 
pulmonary ventilation (VE) at hypoxia from minute 9-10 (h2) (mL/min), breathing 
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frequency (f) during (co) HYPERCAPNIA (breaths/min), breathing frequency (f) during 
(co) HYPOXIA (breaths/min), % change from (co) in frequency (f) under hypercapnia 
conditions to (h2), % change from (co) frequency (f) under hypoxic conditions to (hi), 
breathing frequency (f) at hypercapnia during (hi) (breaths/min), breathing frequency (f) 
5 at hypoxia during (hi) (breaths/min), % change from (co) in frequency (f) under 
hypercapnia conditions to (h2), %change from (co) in frequency (f) under hypoxic 
conditions to (h2), breathing frequency (f) at hypercapnia during (h2) (breaths/min), 
breathing frequency (f) at hypoxia during (h2) (breaths/min), tidal volume (VT) during 
(co) HYPERCAPNIA (ml), tidal volume (VT) during (co) HYPOXIA (mL), % change 

10 from (co) in Tidal volume (VT) under hypercapnia conditions to (hi), % change from 
(co) in tidal volume (VT) under hypoxic conditions to (hi), tidal volume (VT) at 
hypercapnia (hi) (mL), tidal volume (VT)at hypoxia (hi) (mL), % change from (co) in 
tidal volume (VT) under hypercapnia conditions to (h2), % change from (co) in tidal 
volume (VT) under hypoxic conditions to (h2), tidal volume (VT) at hypercapnia (h2) 

1 5 (mL), and tidal volume (VT) at hypoxia (h2) (mL). 

Physiological determinants can include indices of clinical chemistry and 
hematology associated with normoxic and chronically hypoxic conditions in the serum or 
plasma of an organism such as a mammal. Examples include, without limitation, 
amounts of albumin (g/dL), alkaline phosphatase (U/L), alanine transaminase (ALT) 

20 (U/L), anion gap (mmol/L), aspartate transaminase (AST) (U/L), bicarbonate (mmol/L), 
calcium (mg/dL), chloride (mmol/L), cholesterol (mg/dL), creatinine (mg/dL), eosinophil 
(absolute counts in terms of 1000 cells/ (iL), globulin (g/dL), glucose (mg/DI), plasma 
hematocrit (%), hemoglobin (g/dL), lymph (absolute count in terms of 1000 cells/jxL), 
mean corpuscular hemoglobin concentration (pg), mean corpuscular hemoglobin 

25 concentration (g/dL), mean corpuscular volume (fL), and amounts of monocytes 

(absolute count in terms of 1000 cells/^L), phosphorus (mg/dL), platelet count (in terms 
of 1000 cells/uL), potassium (mmol/L), red blood cell (1 X 10 6 /uL), segmented 
neutrophils (in terms of 1000 cells/jiL), sodium (mmol/L), total bilirubin (mg/dL), total 
protein (g/dL), urea nitrogen (mg/dL), and white blood cell count (in terms of 1000 

30 cells/uL). 
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Physiological determinants also can include histological characterization of 
tissues under various physiological conditions, for example, normoxic, hypoxic, or high 
or low salt conditions. Examples include, without limitation, general anatomical 
measurements, measurements derived from medical imaging modalities, or 
quantifications of biomarkers commonly used in disease diagnostics of various tissues, 
for example, those from the aorta, microvasculature, stomach, breast, testes, ovaries, 
bone, lymphocytes, heart, kidney, lung, intestinal, brain, liver, pancreas, and prostate. 

Physiological determinants also can be specific to other disease conditions. For 
example, physiological determinants such as tumor size, cell type, tests of cell type, 
cell/tumor response to various agents, tumor location, primary site of tumor, secondary 
sites of tumors, genes associated with cancer, and microarray patterns of gene expression 
associated with each stage of cancer as well as those described above can be used to 
assess a cancer condition. 

Correlation values and correlation matrices 

As used herein, the term "correlation value" refers to a mathematical relationship 
between two physiological determinants calculated using statistical methods. Standard 
mathematical and statistical methods can be used to determine correlation or other 
statistical or quantitative measures used to characterize relationships between two or 
more physiological determinants. Linear, polynomial, and multiple regression analysis, 
as well as covariance analysis, T-test, and mathematical (linear or non-linear functional 
relationships) are examples of methods that can be used to determine statistical or 
mathematical measures that quantitatively relating two or more determinants. Both 
parametric (model-based) analytical methods (e.g. Pearson correlation coefficient, 
regression methods, mathematical functional relationships) and non-parametric analytical 
methods (e.g. Spearman correlation coefficient, Z-scores, and Wilcoxian rank sum) can 
be used. To obtain a correlation value between two physiological determinants such as 
mean arterial pressure (MAP) and heart rate, for example, the MAP and heart rate for all 
individuals in a study are measured. From the measured values of MAP and heart rate, a 
correlation value, a quantitative measure of the relationship between MAP and heart rate, 
can be determined using the formula: 
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Cxy = S(Xj* Yi) 

(SQRT [Z(Xi*Xj) * 2(Ys*Yi)] ) 

where Xj represents MAP, Yi represents heart rate, and "i" can be 1 to "N". "N" 
5 represents the size of the population being studied. For example, if 100 patients are used 
in the study, then N=100. The 100 values of MAP and 100 values of heart rate are 
presented as 100 pairs of (Xj,Yj). C xy is the correlation value of the MAP and heart rate. 

The correlation value between the two physiological determinants obtained using 
the above formula is the basis of assigning a color to the correlation matrix. The larger 
10 the number of determinants, the larger the correlation matrix, and the more quantification 
required. For M determinants, the matrix is size M, and the number of determinants to be 
calculated is M * M/2. In every case, the mathematical or statistical quantification can be 
ill normalized in such a way as to allow colorization in a consistent manner. All values, for 

J ! f example, are normalized into the range of-1 to 1. This allows for using the same non- 

P 1 5 numerical indications of degrees of correlation. 

j : 1 Other quantifications between any two determinants (Xj, Yj) can be used in 

4* the same manner as the correlation matrix. For example, if the relationship between the 

s 

CI two determinants is non-linear, for example exponential, then the correlation measure 

pi would not provide the most accurate method in characterizing the relationship between X 

□ 20 with Y. Rather, a mathematical model (Y =EXP(B* X) would be the best approach to 
ri characterize the relationship. In this case, the best-fit estimate (based upon a nonlinear 

regression between Y and X) of B would represent the quantified relationship. Hence, 
rather than correlation coefficient, the profile matrix would have the estimate of B in the 
cell representing the relationship between X and Y. 
25 Any combination of correlation coefficients, best-fit estimates, or other 

appropriate measures of the relationships between the variables can be used. For 
convenience, all such measures are referred to herein as "correlation values" between the 
relevant physiological determinants. In some cases, more than one quantification may be 
needed to represent the relationship between two determinants. For example if Y = mX + 
30 b, linear regression would provide the two measures (m, b) representing the quantified 
relationship between Y and X. In this case, more than one "cell" of the profile matrix is 
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required to represent the "profile" attributed to the two determinants (X and Y). 
Whatever the approach to extract quantified measures of the relationships between the 
determinants, in every case there is a simple approach to assigning those measures to the 
correlation matrix, assigning a color or other graphical representation to the measure in 
5 the matrix, and thus creating a colorized or other graphical representation of the 
physiological profile. 

Correlation values can be any value between 1 and -1 inclusive. For example, 
correlation values can be -1, -0.99, -0.9, -0.88, -0.8, -0.77, -0.7, -0.66, -0.6, -0.55, -0.5, - 
0.44, -0.4, -0.33, -0.3, -0.22, -0.2, -0.11, -0.1, 0, 0.1, 0.11, 0.2, 0.22, 0.3, 0.33, 0.4, 0.44, 
10 0.5, 0.55, 0.6, 0.66, 0.7, 0.77, 0.8, 0.88, 0.9, 0.99, 1, or any value in between. 

Once determined, correlation values for all possible pairs of determinants within a 
set of determinants can be presented on a correlation matrix. A "set" of physiological 
J3 determinants is a group of determinants that can be associated with a particular 

m physiological condition. A correlation matrix can be depicted, for example, as a two- 

□ 1 5 dimensional graph in which the determinants are ordered along the X and Y axes (e.g., 

jfl sides, bottom, and top of a two-dimensional array; see, e.g., Figure 3). Correlation values 

^ are placed in locations within the matrix equivalent to locations specified by particular 

□ coordinates (Xs 5 Ys). Determinants can be ordered, i.e. clustered, in a number of non- 
iSI random ways using a clustering method. Determinants can be clustered using known 
Q 20 physiological, biochemical, or functional relationships. For example, all determinants 
U related to a particular biochemical pathway can be clustered next to each other, while all 

determinants related to a biological function, e.g. renal blood flow, can be clustered 
together. As used herein the term "functionally clustered" refers to the ordering of 
determinants based on known physiological, biochemical, or functional relationships. 

25 Determinants also can be clustered using purely statistical or mathematical methods 
involving models (parametric methods) or without models (non-parametric methods). 
Examples include, without limitation, hierarchial, self-organizing maps (SOMs), or 
principal component analysis. Standard statistical methods are described in Everitt, B.S. 
(1993) Cluster Analysis, 3rd Edition, Edward Arnold, Ltd., London, UK; SAS/STAT 

30 User's Guide (1990) Version 6, Fourth Edition, Volume 1, pages 519-614; and 

SAS/STAT User's Guide, 1990, Version 6, Fourth Edition, Vol 2, pages 1614-1631. As 
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used herein, the term "algorithm clustering" refers to clustering determinants using a 
purely mathematical or statistical method. A correlation matrix in which the 
determinants are clustered using functional or statistical methods is herein referred to as a 
"clustered correlation matrix" or a "physiological profile." 

Correlation values can be presented on a clustered correlation matrix as numeric 
values. Correlation values also can be presented in any manner that facilitates visual 
interpretation. For example, a color scheme in which a particular color represents a 
particular degree of correlation can be used. Other types of designations effective in 
differentiating highly negative or positive, moderately negative or positive, or low 
correlation values also can be used, for example shading, stippling, or cross-hatching. 

In contrast to the presentation of correlation values and relationships on a visible 
clustered correlation matrix/physiological profile, the generation of the physiological 
profile can be performed by a computer such that only differences in correlation 
structures under different conditions, or shifts in correlations determined from comparing 
profiles obtained in response to a challenge or over time, are reported to the 
experimenter. In this embodiment, determination of correlation values, statistical 
analyses, phenotypic clustering, and identification of correlation structures or shifts in 
correlations are performed in silico. 

Applications 

Physiological profiling is a method of capturing complex physiological processes 
that contribute to normal and pathological states of an organism. Since physiological 
profiling reveals relationships between physiological determinants, a physiological 
profile generated using determinants related to a complex physiological system can be 
used to capture the efficiency and status of that particular system. 

In one embodiment, physiological profiling can be used to follow development of 
an organism over time. An organism can be subjected to physiological profiling at 
various points in its life cycle. The resulting physiological profiles can be correlated with 
other aspects of the organism's development such as physical, mental, and physiological 
development as well as aging. The resulting physiological profiles also can be correlated 
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with the health status of an organism as well as with susceptibility to infections and 
development of disease conditions. 

In another embodiment, physiological profiling can be performed for large 
populations. Resulting physiological profiles can be used in conjunction with, as 
replacements for, or as supplements to, existing actuarial tables. An individual's profile 
can be used for predicting life expectancy by linking with actuarial tables. 

In another embodiment, physiological profiling can be used as a diagnostic 
method. For example, healthy organisms and those exhibiting, or predisposed to 
developing, a disease condition can be distinguished by physiological profiling based on 
differences in relationships, i.e. correlations, among mechanistically relevant 
physiological determinants. To distinguish a healthy organism from an organism having 
a disease condition, the physiological profiles of organisms or groups of organisms 
representative of normal and disease conditions are determined. The resulting 
physiological profiles representing a normal and a diseased condition can be used for 
diagnostic purposes. 

In another embodiment, physiological profiling can be used to capture 
physiological states of multifactorial diseases. As used herein, the term "multifactorial 
disease" refers to a disease associated with multiple genetic loci as well as environmental 
factors. Examples of multifactorial diseases include, without limitation, obesity, 
hypertension, end stage renal disease, and growth defects. Multifactorial diseases also 
include heart conditions such as myocardial infarction, left ventricular hypertrophy, 
congestive heart failure; diabetes; cancers such as leukemia, lymphoma, and myeloma; 
autoimmune diseases such as lupus, multiple sclerosis, rheumatoid arthritis, type 1 
diabetes mellitus, psoriasis, thyroid diseases, systemic lupus erythematosus, scleroderma, 
celiac disease/gluten sensitivity, and inflammatory bowel diseases; and mental illnesses 
such as schizophrenia, bipolar depression, and Parkinson's disease. 

Typically, the patient population for a particular multifactorial disease, including 
those described above, is heterogeneic in that the clinical presentation of the disease 
condition varies among individuals of the population. Physiological profiling can be 
used to partition heterogeneity, i.e., to reduce the heterogeneous patient population 
exhibiting a multifactorial disease into more homogeneous subclasses of the 
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multifactorial disease. As used herein, the term "homogeneic subclass" refers to a 
subclass of a multifactorial disease population consisting of members whose clinical 
presentation of the disease is more similar to each other than to the clinical presentation 
of the disease in members belonging to another homogeneic subclass of the same 
5 multifactorial disease. 

To identify homogeneic subclasses of a given multifactorial disease, physiological 
profiling of the patient population is performed. Differences in correlation patterns 
identified from physiological profiles can be used to assign patients to homogeneic 
subclasses such that members of each subclass have a physiological profile distinct from 
10 patients in another subclass. Once the homogeneic subclasses of a multifactorial disease 
have been identified, a new patient can be diagnosed as belonging to a particular 
homogeneic subclass by physiological profiling and comparison of the new patient's 
y= ; physiological profile with physiological profiles representative of the different 

^ homogeneic subclasses. The ability to diagnose patients as belonging to particular 

Cj 1 5 homogeneic subclasses of a disease is useful for determining optimal therapeutic 
I ''i regimens. This is because different therapeutic methods can vary in effectiveness 

»> between subclasses. The ability to diagnose a patient as belonging to a particular 

q homogeneic subclasses is also useful for determining prognosis, as particular homogeneic 

fi subclasses may have better survival or other clinical outcomes compared to other 

□ 20 homogeneic subclasses. 

rr In another embodiment, physiological profiling can be used to determine risk 

factors associated with developing a particular disease condition. For example, the 
physiological profiles of patients predisposed to hypertension can be compared to the 
physiological profiles of those not predisposed to hypertension. Correlation patterns 

25 associated with various degrees of predispositions can be identified, and the 

corresponding physiological profiles representative of various risk groups can be used to 
determine the risk group to which a new patient belongs. This is done by comparing the 
new patient's physiological profile with those profiles representing various risk groups. 
In embodiments in which the physiological processes of one organism are 

30 compared to representative physiological profiles in order to, for example, predict 

outcome of a therapy (drug, surgical, biopharmaceutical), determine a prognosis once the 
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disease is identified, or determine an initial prediction of predisposition (actuarial 
assessment of a person's health), a modified method of generating a physiological profile 
is used. In the case of an individual organism, the profile will be produced by assessing 
the relative distance of the physiological determinant value is from the mean population 
value of the same physiological determinant. This distance will then be used in a manner 
similar to the correlation value. The overall pattern of the profile then is analogous to the 
correlation matrix. Prognosis, diagnosis and predisposition can then be determined 
empirically by the similarity or difference in the individual's profile versus other patients' 
known outcomes with similar profiles or the population average. This predictive nature 
of the profile can be used for various organisms, for example, humans. 

Physiological profiling can be used in combination with genetic linkage analysis 
to identify loci associated with different clinical presentations, i.e. symptoms or 
manifestations, of a multifactorial disease. Populations of patients having a multifactorial 
disease condition typically exhibit heterogeneous clinical presentations. Physiological 
profiling allows the heterogeneous patient population to be partitioned into more 
homogeneous subclasses. Genetic linkage analysis can identify chromosomal regions 
that are associated with particular phenotypic presentations. Combining physiological 
profiling with genetic linkage analysis data allows for identification of multiple genetic 
loci that may give rise to similar clinical presentations. 

In another embodiment, physiological profiling can be used as a comprehensive 
approach to characterizing the influences of particular genomic regions on the 
relationships among pathways within complex physiological processes. Genetic linkage 
analysis alone reveals the direct influences of genes on the mechanisms measured by the 
mapped phenotypes. The influences of genes on mechanisms measured by the mapped 
phenotypes represent first order linkage. Physiological profiling allows for identification 
mechanistic relationships among pathways associated with complex physiological 
processes. When genetic linkage analysis is combined with physiological profiling, the 
effects of genotype on relationships among pathways within complex physiological 
process can be determined. Thus, combining genetic linkage analysis with physiological 
profiling provides a means to relate genetic information with functional pathways. 

Physiological profiling also can be combined with expression profiling, either 
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alone or in combination with genetic linkage analysis, to perform functional genomics. 
Expression profiling, as described, for example in U.S patent numbers 6,251,601; 
5,800,992; and 5,445,934 can be used to identify genes that are expressed under 
particular conditions. Genetic linkage analysis identifies locations of the genome that are 
associated with particular phenotypic determinants. Physiological profiling identifies 
relationships among phenotypic determinants. Knowledge of the expression profiles of 
individual genes, their locations on chromosomes, and their effects on relationships 
among functional pathways within complex physiological processes can provide 
profound insights into the biology of organisms. 

The invention will be further described in the following examples, which do not 
limit the scope of the invention described in the claims. 

EXAMPLES 

Example 1 - Animals used in a study on the genetic basis of hypertension 
F2 progeny rats derived from an intercross of an inbred hypertensive rat and a 
normotensive rat were used. The inbred hypertensive rat was a Dahl salt sensitive rat 
(SS/JrHsdMcw), and the inbred normotensive rat was a Brown Norway rat 
(BN/SsNHsdMcw). Two hundred and twelve F2 rats (113 males and 99 females) were 
extensively phenotyped for 239 mechanistically relevant cardiovascular, neuroendocrine, 
and renal phenotypes, including a number of cardiovascular stressors, both dietary and 
pharmacological, as described in Examples 2, 3, 4, and 5. 

Example 2 - Phenotyping protocol for conscious animals at high and low salt intakes 
Rats were maintained on a high salt diet (8 % salt) from the age of 9 to 13 weeks. 
During the fourth week of the high salt diet, arterial pressures of un-anesthetized rats 
were measured for three hours each day for three days. All blood pressure (BP) 
measurements were made with the animals unrestrained in their home cages as described 
previously (Cowley Jr. et al (2000) Physiological Genomics 2:107-1 15). Implanted 
arterial catheters were used in determining arterial pressures. Data were collected at a 
rate of 1 00 Hz and reduced to one-minute averages; data for time series analysis were 
reduced to one-second averages. At the end of the third high salt day, animals were salt- 
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depleted and placed on a low salt diet. One and a half days following furosemide- 
induced salt depletion and switching to a low salt diet, arterial pressure responses were 
determined. The day-night light cycle for all rats ran from 2:00 AM (lights on) to 2:00 
PM (lights off) throughout the study. 

Blood pressure data for high salt day 1 (BP1) consisted of baseline measurements 
of heart rate and systolic, diastolic, and mean arterial pressures measured from 9:00 AM 
to noon. 

Blood pressure data for high salt day 2 (BP2) consisted of measurements of heart 
rate and systolic, diastolic, and mean arterial pressures obtained for the inactive (lights 
on) and active (lights off) phase. Data for the inactive phase (baseline data) were 
obtained from 9:00 AM to noon as was done for high salt day 1 . Data for the active 
phase were obtained from 2:00 PM- 6:00 PM. All blood pressure data on this day were 
collected for time-series analysis. A 24-hour urine collection was started in which urine 
volume as well as sodium, potassium, protein, and creatinine levels were determined. 

Blood pressure data for high salt day 3 (BP3) consisted of baseline measurements 
of heart rate and systolic, diastolic, and mean arterial pressures measured from 9:00 AM 
to noon. Following the baseline measurements, a blood sample (500 |iL) was drawn for 
determination of plasma renin activity and creatinine, plasma protein, and hematocrit 
levels. Following the blood draw, an injection of furosemide (10 mg/kg) was given 
intraperitoneally (ip) to salt deplete the animals. Following the furosemide 
administration, the animals were switched to a low salt diet (0.4 % salt). 

Blood pressure data for salt-depleted-day 4 (BP4) consisted of measurements of 
heart rate and systolic, diastolic, and mean arterial pressures measured from 9:00 AM to 
noon in the salt depleted state. These measurements were followed by a stress test. The 
stress test consisted of delivering two alerting stimuli five minutes apart; each alerting 
stimulus was 2 milliamps for 0.3 seconds. The change in mean arterial pressure, the time 
to peak, and the time to 90 % recovery in response to the stress test were determined. 

Blood pressure data for salt depleted-day 5 (BPS) consisted of measurements of 
heart rate and systolic, diastolic, and mean arterial pressures determined from 9:00 AM to 
noon in the salt depleted state. Following the recording period, a 1.0 mL blood sample 
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was taken for determination of plasma renin activity; white blood cell count; and 
triglycerides, total cholesterol, HDL, creatinine, and hematocrit levels. 

Example 3 - Phenotyping protocol for renal and peripheral vascular reactivity in 
5 anesthetized animals 

Rats were anesthetized with 30 mg/kg of ketamine and with 50 mg/kg of Inactin 
administered intraperitoneally. Catheters were implanted in the femoral artery and vein, 
and an electromagnetic flow probe was placed on the left renal artery via a midline 
incision. An intravenous (iv) infusion (50 (iL/min) of isotonic saline containing 1 % 
1 0 bovine serum albumin was performed to replaced fluid loss. After a 45-minute 

equilibration period, control values of arterial blood pressure and renal blood flow (RBF) 
C5 were measured for 15 minutes. Next, animals were given iv infusions of angiotensin II 

(20, 100, 200 ng/kg/min) and norepinephrine (0.5, 1, 3 (ig/kg/min) for 5 minutes after 
JLJ which renal and peripheral vascular responses were determined. Following recovery of 

fU 15 pressure to baseline values, animals were given two successive doses of acetylcholine 
]= (ACh) (0.1 and 0.2 jag/kg/min as bolus doses) after which renal vascular and systemic 

^ arterial responses were measured. To determine the contribution of nitric oxide to basal 

y3 renal vascular tone, 5 mg/kg of nitro-L-arginine methyl ester (L-NAME) was 

?Z administered as an iv bolus, and then renal blood flow and renal resistance were 

W 20 determined. After 10 minutes of equilibration, the degree of blockade of the synthesis of 
nitric oxide produced by L-NAME was determined by administration of a repeat infusion 
of the same two doses of Ach, and renal blood flow and renal resistance were examined. 



Example 4 - Collection of tissue samples for morphometric measurements and histology 
25 To assess the degree of cardiac and renal hypertrophy, heart and kidneys were 

removed, stripped of surrounding tissue, and weighed using a digital top loading 
Sartorius balance. For histological analysis, the right kidney was fixed by immersion in 
1 0 % buffered formalin, embedded in paraffin, and the prepared sections were stained 
with hematoxylin and eosin as well as Periodic acid Schiff (PAS). These sections were 
30 evaluated for mean glomerular diameter and the degree of focal glomerulosclerosis. The 
degree of focal glomerulosclerosis was used as an index of glomerular injury. 
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Example 5 - Development of a genetic linkage map determinant phenotypes of blood 

pressure 

To obtain a comprehensive picture of the genomic regions that are linked to blood 
pressure determinants, a genetic linkage map of measured and derived determinant 
phenotypes obtained as described in Examples 2, 3, and 4 was generated. These 
phenotypes represented critical elements of neuroendocrine, vascular, and renal 
functions. Table 1 below summarizes the physiological determinants used this study. 
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Determinant phenotypes were tested for normalcy as described in Cowley Jr. et 
al (2000) Physiological Genomics 2: 1 07- 1 1 5. The 1 66 phenotypes that passed a test for 
being distributed as a normal random variable were analyzed via parametric methods in a 
genome scan using MAPMAKER/QTL (as described in J. P. Rapp (2000) Physiol Rev. 
5 80:135-172 and Hollenberg et al (1978) Medicine 57:167-178). The remaining 73 
phenotypes that did not fulfill the requirements for parametric analysis were analyzed 
using a non-parametric mapping algorithm (MAPMAKER/QTL version 1.9b). Distances 
between loci were calculated based on the Haldane algorithm. An average spacing of 
markers of 1 0 cM was used. 
10 To determine the threshold for suggestive and significant linkage, a permutation 

test was performed. The permutation test consisted of 5000 random assignments of 
genotypes with phenotypes. These results confirmed that the LOD thresholds of 2.8 and 
4.3 for suggestive and significant linkages, respectively, set by Lander and Kruglyak 
jjf (Loscalzo et al (1995) Progress in Cardiovascular Diseases 38:87-104) were 

P 1 5 appropriate for this study despite the large number of phenotypes tested. 

Eighty-one phenotypes had either a parametric LOD score > 2.8 or non- 
parametric LOD score > 3.5; 18 parametric phenotypes had LOD scores between 2.5-2.8; 
and 26 phenotypes were functionally related to blood pressure. From these 81 parametric 
and non-parametric phenotypes, 96 QTL were identified of which 69 had an LOD score 
20 of > 2.8 and 25 had a LOD score of > 3.5. The 96 QTL identified in the autosomal 

genome of 1 13 male progeny from an SS/JrHsd/Mcw x BN/SsNHsd/Mcw intercross are 
shown in the genetic linkage map of Figure 1. 



■S3ET 
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Example 6 - Results of genetic linkage analysis 
25 In general, QTL for blood pressure were clustered in discrete regions on rat 

chromosomes 1, 2, 3, 7, and 18. These clusters consisted of six or more QTL with 
overlapping 95 % confidence intervals. In four of the five clusters, the determinant 
phenotypes were independent, indicating that these four clusters represented separate 
genes rather than a pleiotropic effect. In the fifth cluster, on chromosome 1 8, significant 
30 correlations were found among the determinant phenotypes. These phenotypes could be 
divided into three functional groups that include phenotypes associated with: vascular 
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reactivity, plasma lipid concentration, and renal function. The clustering of QTL 
associated with phenotypes belonging to distinct functional groups suggests the presence 
of a functional cassette as has been observed for QTL in agriculture and biomedical 
research. (See Thumma BR et al (2001) J. Exp. Bot Feb, 52, 203-214; Miner LL, M. RJ 
(1995) Psychopharmacology (Berl) 1 17, 62 - 66; Wakeland et al (1997) J Clin Immunol 
17, 272-281; and Nadeau et al (2000) Nat Genet 25, 381-384). 

More specifically, QTL for MAP and RBF responses to ACh were found on 
chromosome 10 (DIOMghH); this contributes 17 % to the variance of the pressure and 
RBF in the F2 population. When the contribution of DIOMghH to genetic variance was 
removed using the fix command in MAPMAKER (Lander et al (1987) Genomics 1:174- 
181), additional loci on chromosomes 4 (D4Mit2) and 12 (D12Mit7) were found to 
contribute to the ACh response. The loci on chromosome 10, 4, and 12 were known to 
harbor genes for nitric oxide synthesis, i.e. NOSIII on chromosome 4, NOSH on 
chromosome 10, and NOSI on chromosome 12. While increases in the synthesis of nitric 
oxide mediate much of the vasodilator response to ACh (Loscalzo et al (1995) Progress 
in Cardiovascular Diseases 38:87-104), until now, the vasodilator response has not been 
shown to be associated with all three nitric oxide synthases. 

Example 7 - Development of physiological profiling for studying physiological responses 

Phenotyping and genetic mapping data obtained as described in Examples 2-6 
were used in developing a new analytical strategy - physiological profiling. Genetic 
mapping of determinant phenotypes associated with hypertension resulted in 
identification of QTL that, in many cases, were clustered in the same region of the 
genome. To understand the chromosomal clustering of the phenotypes, a correlation 
matrix was constructed. The correlation matrix consisted of correlation values 
determined by linear regression analyses for all pairs of phenotypes. Each correlation 
value reflected the relationship between two phenotypes. For ease of visual analysis, 
correlation values ranging from 1 to -1 were presented on the matrix using a color 
scheme. Figure 2, for example, is a colorized correlation matrix of the genetically 
mapped phenotypes of BN rats. Phenotypes were organized in a random order along the 
top and side of the colorized correlation matrix. To capture the complex interactions 
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among phenotypes, a second analytical procedure was performed to cluster phenotypes 
into meaningful groups. The clustering procedure was performed using known functional 
or physiological relationships (functional clustering). In functional clustering, for 
example, related phenotypes such as RBF responses to agonists were placed next to each 
other. Figure 3 is a functionally clustered correlation matrix, i.e. a physiological profile, 
consisting of the same phenotypes used in Figure 2. In the physiological profile depicted 
in Figure 3, the phenotypes were clustered based on Guyton's model of blood pressure 
control (Guyton, A.C. (1972) Monograph). 

To validate the use of functional clustering, a physiological profile in which 
phenotypes were clustered using a two-step clustering algorithm that was independent of 
function was generated. (See, Everitt, B.S. (1993) Cluster Analysis, 3rd Edition, Edward 
Arnold, Ltd., London, UK; SAS/STAT User's Guide (1990) Version 6, Fourth Edition, 
Volume 1, pages 519-614; and SAS/STAT User's Guide, 1990, Version 6, Fourth 
Edition, Vol 2, pages 1614-1631.) The physiological profile generated using the two-step 
algorithm clustering method was compared with one generated using functional 
clustering based on Guyton's model of blood pressure control (Figure 3). The two 
matrices were overlayed to form a composite profile (Figure 4). Comparison of the 
composite profile of Figure 4 with the physiological profile generated using Guyton's 
model of blood pressure control (Figure 3) showed that they were very similar. These 
results suggest that a correlation matrix generated by functional clustering is useful for 
establishing a "profile" of physiological function. 

Example 8 - Comparison of the physiological profiles of the parental BN and SS rats 
Correlations among phenotypes of the parental BN and SS rats that corresponded 
to the mapped phenotypes of the F2 intercross were analyzed by physiological profiling. 
The physiological profiles, shown in Figure 5, were generated from 50 parental BN rats 
(top right triangle) and 50 parental SS rats (bottom left triangle). Phenotypes were 
functionally clustered, i.e. using knowledge of physiological function and without the aid 
of the correlation matrix. Comparison of the physiological profiles of these two strains 
shows clear differences in a number of correlations. For example, low positive 
correlations were observed in the RBF phenotypes of the BN rats, while significant 



43 



Docket 13482-< 




negative correlations among the same phenotypes were observed in the SS rats. These 
results are consistent with the finding that salt sensitive SS rats cannot regulate renal 
resistance and blood flow as well as BN rats in response to renal perfusion pressure. The 
agreement with prior observations illustrates the effectiveness of physiological profiling. 
Second, the two prominent clusters outlined in red in both strains represent results of 
dose response curves of acetylcholine, angiotensin II, and norepinephrine. These clusters 
indicate that both strains have similar physiological responses to these pharmacological 
agents. 

Example 9 - Comparison of the physiological profiles of BN rats and all F2 intercross 
progeny rats generated by functional and algorithm clustering 
The physiological profile of BN rats was compared with the physiological profile 
of all F2 intercross progeny (see Figure 6). Two physiological profiles of the parental 
BN rats, one generated using functional clustering and the second using a clustering 
algorithm, were compared by overlaying. Similarly, two physiological profiles of the F2 
intercross progeny, one generated using functional clustering and the second using the 
same clustering algorithm as in Example 7, were compared by overlaying. The resulting 
composite profiles of the BN rats (top right and above diagonal) and the F2 progeny rats 
(bottom left and below the diagonal) are shown in Figure 6. A blending of profiles in the 
F2 intercross progeny was observed when the composite BN and F2 progeny profiles 
(Figure 6) and the SS profile in Figure 5 were compared. Although functional clustering 
and clustering by purely statistical methods yielded similar results as indicated by the 
composite profile, the physiological profile generated by statistical methods revealed two 
clusters of traits (clusters 12 and 14) that were not known before. 

Example 10- Effects of allelic substitution on the physiological profiles 
Physiological profiling was used to assess how a group of F2 animals respond to a 
salt challenge. Figures 7 A and 7B are comparisons of the physiological profiles of the 
entire F2 population with the F2 animals that have QTL that protect against a salt load. 
The F2 animals that were protected against a salt load were those that fell into the 10 % 
tails of a distribution after a salt challenge. In general, low correlations between 
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phenotypes were observed in the physiological profile of the entire F2 population. In 
contrast, strong positive correlations between some phenotypes were observed in the 
physiological profile of animals that were protected against a salt challenge. Therefore, 
physiological profiling is effective in capturing differences in physiology between 
5 distinct groups. 



Example 11 - Physiological profiling of blood pressure determinant phenotypes 
Since genetic mapping results of Example 6 demonstrated that the vasodilator 
response to Ach is associated with loci containing three nitric oxide synthases, the impact 
10 of BN and SS alleles of all three NOS genes on the mapped phenotypes was examined by 
physiological profiling. This also allowed for assessing the systems biology of the other 
mapped cardiovascular phenotypes. 

CJ 

J3 The mapped phenotypes were analyzed by physiological profiling using 

~iz functional clustering based on Guyton's model of blood pressure control. Figure 8A is 

C3 1 5 the physiological profile for F2 male rats that were homozygous SS for D10Mghl4 (the 

\ i J 

i'i flanking marker for NOSH), and Figure 8B is the physiological profile for F2 male rats 

** a that were homozygous BN for D10Mghl4. The correlation patterns were found to be 

f a quite different when the SS and BN profiles were compared. In F2 rats homozygous for 

the SS allele at D10Mghl4 (NOSII), positive correlations were observed among blood 
£3 20 pressures determined immediately before, during, and after the short-term intravenous 
administration of norepinephrine (NE), angiotensin II (Ang II), and acetylcholine (see 
Figure 8A, cells #85-94). In contrast, F2 rats homozygous for the BN allele at 
D10Mghl4 (NOSII) exhibited weak correlations among the same phenotypes (see Figure 
8B). Furthermore, the relationships of measured differences in systolic, diastolic, and 
25 mean arterial pressures before and following infusions of NE (cells #88-92) also were 
different. 

The relationship between MAP before and after infusion of NE in F2 rats carrying 
the BN or SS allele at DIOMghH (NOSII) is further illustrated in Figures 9A and 9B. 
Figure 9A is a graph demonstrating the correlation between MAP before and after 
30 infusion of NE in F2 rats carrying the BN allele (closed circles) and in F2 rats carrying 
the SS allele (open circle) at DIOMghH (NOSII). Although the F2 rats carrying the BN 
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allele (closed circles) at D10Mghl4 exhibited a lack of correlation between arterial 
pressure levels before and following intravenous infusions of NE, a significant positive 
correlation was observed in F2 rats homozygous for the SS allele (open circles) at NOSII. 
Figure 9B is a bar graph summarizing the average levels of MAP before (solid bars) and 
following completion (open bars) of intravenous infusions of three doses of 
norepinephrine in male rats carrying the SS or BN allele at NOSII. Arterial pressures of 
male rats carrying the SS allele at NOSII returned precisely to control levels, while the 
arterial pressures of rats carrying the BN allele fell significantly below control and 
remained at hypotensive levels for as long as 10 minutes. Although it has been reported 
that NO plays a role in the control of vascular tone and blood pressure, the finding that 
NOS is involved in vasoconstrictor agents-induced hypotension has not been reported. 

The physiological profile of the heterozygote at D10Mghl4 exhibited a 
correlation pattern that was intermediate between SS and BN, although closer to BN (data 
not shown). 

Physiological profiles for rats partitioned by alleles on chromosome 4 at D4Mit2 
(NOSIII) demonstrated similar relationships for these traits (#85-94) for SS versus BN 
alleles (see http://brc.mcw.edu/phvprf) . The similarity in physiological profiles of rats 
partitioned by at NOSIII and those of rats partitioned at NOSII indicate similar gene 
effects on overall physiology. 

Example 12 — Other novel relationships derived from physiological profiling 
Positive correlations were observed between the urinary excretion of protein and 
plasma lipid concentration (Figures 8 A and 8B, cells 1 15 vs. 34), and between kidney 
weight and plasma lipid concentrations (Figures 8 A and 8B, cells 1 15 vs. 36) in F2 rats 
homozygous BN for the NOSII allele. Although hyperlipidemia, proteinuria, and renal 
hypertrophy are related symptoms seen in hypertensive and diabetic nephropathy and end 
stage renal disease, the influence of a NOSII genotype on the relationships between these 
indices of renal end organ damage and hyperlipidemia has not been described previously. 
Thus, the present finding that the severity of proteinuria, renal hypertrophy, and 
hyperlipidemia in an F2 population of rats is influenced by the allelic variations of 
inducible NOS (NOSII) gene is novel and creates new directions for further research. 
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Another relationship determined by physiological profiling that was not detected 
by linkage analysis was the strong positive correlation found between Angiotensin II- 
induced reductions of RBF and the chronic urinary excretion of protein in rats 
homozygous BN at the D13Mghl8 allele. This finding contrasts with the significant 
negative correlation observed in F2 homozygous SS rats at this same allele. (Urine 
protein excretion was determined during steady-state conditions of high salt intake, while 
renal blood responses to Angll were determined following a day of salt depletion in 
anesthetized F2 rats using three doses of AngIL) These results provide a genetic basis for 
results of many clinical studies of essential hypertension in which patients have been 
stratified based on their renal vascular responses to Angll as "modulators" or 
"nonmodulators" (Hollenberg et al (1978) Medicine 57:167-178). In 40-50 percent of 
the essential hypertensive population, adrenal and renal vascular responses to Angll are 
not modified by changes in sodium intake as would be expected. These individuals have 
been called "non-modulators" (Hollenberg et al (1978) Medicine 57:167-178). It has 
been documented that non-modulators exhibit a higher percentage of one or both parents 
with hypertension suggesting that this renal abnormality is inherited and linked to the 
development of hypertension (see Hollenberg et al. (1978) Medicine 57:167-178). In the 
present study, SS fed a high salt diet exhibited substantial proteinuria compared to BN 
parental rats. The physiological profile revealed that in F2 rats, it was possible to predict 
the renal blood flow Angll sensitivity based on genotype and protein excretion levels. A 
narrow region on rat chromosome 13 near D13Mghl8 enables a prediction of modulators 
and non-modulators using individual genotype. The only identified gene that maps very 
closely to D13Mghl8 is renin, an obvious candidate for these responses. The present 
result provides a genetic basis for modulators and nonmodulators that could be explored 
further in a genetic rat model of hypertension. 
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Example 13 - Physiological profiles of African American and French Canadian patients 



Physiological profiling was used to assess correlations between phenotypes 
associated with blood pressure in resting and stressed patients with hypertension. 
5 Patients were African American and French Canadian sibling pairs with hypertension. 
Patients underwent an extensive 2-day in-house protocol (see Kotchen et al. (2000) 
Hypertension 36:7-13 and Pausova et al. (2001) Hypertension 38:41-47) at the Medical 
College of Wisconsin or Centre de Recherche, Centre Hospitalier de l'Universite de 
Montreal (CHUM), Montreal, Canada. Figure 10 depicts the physiological profiles 
1 0 generated for the French Canadian and the African American patient populations. A 
comparison of the physiological profiles of the two patient cohorts revealed that the 
sibling-sibling profiles were more similar than the matrices generated when only one of 
the siblings was used. These data indicate that physiological profiling is useful for 



complex physiological interactions into a single image, (2) capture graphically in a single 
image the numerous differences in the cardiovascular system of different rats strains, and 
(3) identify physiological characteristics not evident by genetic mapping data. 



It is to be understood that while the invention has been described in conjunction 
with the detailed description thereof, the foregoing description is intended to illustrate 
and not limit the scope of the invention, which is defined by the scope of the appended 
claims. Other aspects, advantages, and modifications are within the scope of the 
25 following claims. 



with hypertension 




comparing heritable traits. 

Therefore, physiological profiling is a powerful means to (1) summarize the 
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